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Abstract 

We analyze the decay B — ► <pK within the framework of QCD-improved factorization. We found 
that although the twist-3 kaon distribution amplitude dominates the spectator interactions, it will 
suppress the decay rates slightly. The weak annihilation diagrams induced by (S — P)(S + P) 
penguin operators, which are formally power-suppressed by order (AQCD/^b) 2 ) are chirally and 
logarithmically enhanced. Therefore, these annihilation contributions are not subject to helicity 
suppression and can be sizable. The predicted branching ratio of B~ — ► <fiK~ is (3.8 ±0.6) x 10 -6 in 
the absence of annihilation contributions and it becomes (4.3^ 4) x 10~ 6 when annihilation effects 
are taken into account. The prediction is consistent with CLEO and BaBar data but smaller than 
the BELLE result. 
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I. INTRODUCTION 



Previously CLEO has put an upper limit on the decay mode B — > <pK [T|]: 

B(B ± -> 0if+) < 5.9 x HT 6 . (1.1) 
However, both CLEO 0, BELLE || and BaBar have recently reported the results: 

f (5.5i?j±0.6) x 10- 6 CLEO, 
B{B ± -> 07^+) = | (7.7j£! ± 0.8) x lO" 6 BaBar, (1.2) 
I (13.9±i:I±i;|) x 10- 6 BELLE, 

and 

f (5.4±|;? ± 0.7) x 10~ 6 < 12.3 x 10~ 6 CLEO, 
B(B° -> 0K°) = | (8.l+3;i ± g) x 10 _ 6 BaBar; (i.3) 

I < 16.0 x 10- 6 BELLE. 

It is known that the neutral mode B — > <f)K° is a pure penguin process, while the 
charged mode (f)K~ receives an additional (though very small) contribution from the tree 
diagram. The predicted branching ratio is very sensitive to the nonfactorizable effects which 
are sometimes parameterized in terms of the effective number of colors N% \ it falls into a 
broad range (13 ~ 0.4) x 10 -6 for N^ s = 2 ~ oo ||. Therefore, a theory calculation of the 
nonfactorizable corrections is urgently needed in order to have a reliable prediction which 
can be used to compare with experiment. 

A calculation of B — > <pK within the framework of QCD-improved factorization has been 
carried out recently in 0. However, the analysis of is limited to the leading order in l/m& 
and hence the potentially important annihilation contributions which are power-suppressed 
in the heavy quark limit are not included. 

In the present paper we will analyze the decay B — > <pK within the framework of QCD- 
improved factorization. We will study the important twist-3 effects on spectator interactions 
and also focus on the annihilation diagrams which are customarily assumed to be negligible 
based on the helicity suppression argument. However, weak annihilations induced by the 
(S — P)(S + P) penguin operators are no longer subject to helicity suppression and hence 
can be sizable. This is indeed what we found in this work. 



II. GENERALIZED FACTORIZATION 



The effective Hamiltonian relevant for B — > <pK has the form 



H cS (AB = 1) = ^=lv ub v: s 



ci(/i)Oi(/i) + c 2 (/i)0 2 (/i) 
VtbV* ( £ c z (fi)Oi(fi) + c g (fx)O g (fx) 



10 



+ h.c, 



(2.1) 



\i=3 



2 



where 



O x = (ub) v _ A (su) v _ A , 

03(5) = Wv-iE^v-iM. 



^Mv-aZX^Wv-^v-m)' ( 2 - 2 ) 



v+a(v~a)' 



°7(9) = 2( s6 )^E e l'(«'?'i 

9 = ^m 6 ^G^y(l+ 75 )6 



2 = (« a &/3)v„ A (s /3 U a ) v _ A , 
4 (6) = 1 

3 

°8(io) = ^{s a bp) v _ A Y,^'{q'f3q' a , 



V+A 



iv~A)l 



with {qiq 2 ) v±A = 917^(1 ± 75)92, 3 -0 6 being the QCD penguin operators, 7 -O w the 
electroweak penguin operators, and O g the chromomagnetic dipole operator. 

In the generalized factorization approach for hadronic weak decays, the decay amplitudes 
oiB^(f)K read (in units of G F /V2) 



A(B~ K~ 



a 3 + a 4 + a 5 - -(a 7 + a 9 + a 10 ) 



+ a 4 + a w — 2(a 6 + a 8 



B~ 



(m s + m u )(m b + m u ) 



x (B-K-,<f,) 

X {B-,4>K~) 



+V ub V: s a l X^ B -^ 

-v tb v* 



(2-3) 



a 3 + a 4 + a 5 — -(a 7 + a 9 + ai ) 



A 



1 



ml 



^4 - -aio - (2a 6 -as)/ , w 

2 (m s + m d )(m b + m d ) 



x (B u ,<piC) 



where the factorized terms 



= (0|(« S ) v _JO)(ir|(g6),_j5) = 2Um 4) F 1 BK {ml){e* -p B ), 
X (**+) ee &K\(sq) v _ A \0)(0\(qb) v _ A \B) = 2/ B mXW)(^ • P fl ), 



(2.4) 



can be expressed in terms of the form factors Ff K , (for the definition of form factors, 
see ||) and the decay constants f<p and and the nonfactorized contributions parametrized 
in terms of Xi are lumped into the effective number of colors N^ s : 



NfJt 



(2.5) 



so that the effective parameters dj appearing in Eq. (|2.3| ) read 



a 2i — C-2i + 



"C2i-1, 



a 2i-l — c 2i-l + 



-C2i- 



(2.6) 



(Nf) 2i — 11 ^ ^ (iVf ) 2i _ a 

It is known that the parameters a 3 and a 5 depend strongly on N* G , while a 4 is iV^-stable 
(see, for example, |7]|§). Therefore, the prediction of 5 — ► rates is sensitive to N^ s and 
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hence to the nonfactorizable terms Xu it varies from 13 x 1(T 6 to 0.4 x 1(T 6 for Nf ranging 
from 2 to oo [Q. 

Owing to the unknown form factor Aq K \m 2 B ) at large q 2 , it is conventional to neglect the 
annihilation contribution based on the argument of helicity suppression, which amounts 
to having a vanishing form factor Aq K (i7i 2 b ). However, this argument is valid only for 
(V — A)(V — A) interactions but not for (S — P)(S + P) ones. This explains the large 
enhancement factor of m 2 B /(mbm s ) for the penguin contributions [see Eq. ( |2.3j )l- Therefore, 
it is conceivable that the annihilation contribution could be sizable and significant. 



III. NONFACTORIZBALE EFFECTS IN PENGUIN AMPLITUDES 



We next proceed to compute the nonfactorizable effects in the QCD-improved factoriza- 
tion approach. For simplicity we will neglect the light quark masses. In the chiral limit, 
kaon is massless, but the meson has a finite mass. We consider the vertex corrections and 
hard spectator interactions depicted in Fig. 1 as well as the annihilation diagrams shown 
in Fig. 2. Recently we have analyzed B — > J/ipK decays within the framework of QCD 
factorization |10|| . The study of B — > <j)K is quite similar to the J/ipK mode except for the 
absence of weak annihilations in the latter. The reader is referred to [ll| for details. The 
resultant amplitudes are 



A{B' 



K~ 



a 3 + a 4 + a 5 - -jifli + a 9 + a 10 ) 



x (B-K-,<t>) 



A(B° 



K 



+ 



(c 3 + cv)Al f + (c fi + c 7 )A 2 nf + (c 6 + c 8 + -(c 5 + c 7 ))A f 



a 3 + a 4 + a 5 — -(a 7 + a 9 + a 10j 



X (B°K°, 



+ 



V ub V: s c 2 A l nf1 
(3.1) 



(c 3 - )fv)Al f + (cs - ^c 7 )A 2 nf + (c 6 + ^c 5 - ^c 8 - ^c 7 )A f 



where 



C4 a s Cf 

a3 = C3 + iv; + 47iv: C4 
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c 3 a s C F \ 



14 ) rrib 

f^)-121nA + /; + / „ 
V 14 I m b 



\ x L G{s u ) + -^G{s c ) J + Yl ( c 4 + c 6 + \e q {c s + c w ))G{s q ) +c g G g \, 

V A * A t I q=u,d,s,c,b 1 ) 



c 6 a s C F 



a 7 = c 7 + 



c 8 



a., 



iV„ 4vr JV, 



■C8 



_ u\n-^ + f I + f n 

m b 

- 121n-^- + f! + f n 



(3.2) 



a 

9^ 



N c C e , 



4 





(6) 





(d) 





(«0 



(/) 



FIG. 1. Vertex and spectator corrections to — ► <f>K . 



o 9 



aio 



. cio a s C F 

C9+ A^ + 4^iv: Ci0 

c 9 a s CV 

Cl ° + iv: + ^ iv: C9 
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- 12 In 


JL 


14; 




m b 


ISA 


- 12 In 


jj_ 


14 J 




m b 



—N c C e 
9tt 



+ // + fn 



a: 

9tt' 



a. 



In Eq. ( p.2|) , the upper entry of the matrix is evaluated in the naive dimensional regulariza- 
tion (NDR) scheme for 75 and the lower entry in the 't Hooft-Veltman (HV) renormalization 
scheme, Cf = {N% — 1)/(2N C ), s q = m/i/ml, X q i = V q 'bV*, s , and a the electromagnetic 



q'o • q'si 

are 



fine-structure coupling constant. The other terms in (|3J 

G(s) = - - - In + 4 f 1 d£ ) / duuil - u) ln[s - u(l - u)(l - 0], 
' ' Jo Jo 



2 4 n At 
In — + 4 

3 3 m& 



(c 2 + 



Cl_ 



(3.3) 



where 3>^(£) is the light-cone distribution amplitude (LCDA) of the <fi meson, which will be 
discussed shortly. B — > <pK do not receive factorizable contributions from a§ and ag except 
for the annihilation topologies. The nonfactorizable annihilation contributions A^ and the 
factorizable annihilation amplitude Af will also be elucidated on below. 

Note that the effective parameters d{ appearing in Eqs. ( |3.2|) are renormalization scale 
and 75-scheme independent. Since only one gluon exchange has been considered in the 
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annihilation diagrams (see Fig. 2), one may wonder the scale and scheme dependence of the 
annihilation amplitude given in Eq. Fortunately, as we shall see below, the annihilation 

contribution is predominated by the penguin effects characterized by the parameters = 
cq + |c 5 and a 8 = c 8 + |c 7 multiplied by fi x . It turns out that the scale dependence of a 6 
and as is canceled by the corresponding dependence in /i x owing to the runing quark masses. 
Consequently, the annihilation amplitude is essentially scale independent. 
The hard scattering kernel // appearing in Eq. ( |3.2| ) reads 



fi 



'3(1-20 



lnf 



+ 



(-Lzl 



Sm H : h 3 m(l — z) 



z 2 i In zi + 



z 2 e 2 [ln(l -z) -in] 
[l-z(l-0? 



(3.4) 



where z = m^/m B . For completeness, we have included the mass corrections to //,[] 
though such corrections are very small. In the — > limit, // has the same expression as 
that in B — > tttt decay [11]], as it should be. The hard scattering kernel fu arises from the 
hard spectator diagrams Figs, le and If and has the form [ TO | 

a a (p h ) 4-7T 2 f K fB 



a s (/x) N c F* K {m 2 



13 



1 — z Jo p JO C, 



X 



1 drj 

o fj 



where 



2 /*x(^/») 
m B 



2m 2 K 



l-z] 



Qrj _ 



(3.5) 



-Mm) 
fi 



m s {p) + m u {p) JK 
is proportional to the quark condensate, the B meson wave function $f is defined by |T 

(o\U^b(o)\B( P ))\ x+=x±=0 ^ 

with n_ 



(3.6) 



[{jt> + m B )^] Pl / dpe-^-[^(p)+^(p)) 7a , (3.7) 



'1, 0, 0, —1), and $^ is a twist-3 kaon LCDA defined in the tensor matrix element 



(K-{P)\s{0)^ 5 u{x)\0) 



fKm 2 K 



m,. 



1 



6 m s 



m s + m u 
m K 

dfje if > Px <S>*{fi). 



(3- 



*Eq. (|J) can be obtained from Eq. (19) of g2| or from Eqs. (2.22) and (2.23) of JTrJ by 
neglecting the £ 2 terms arising from the transverse wave function and applying the relation 
FQ K (m 2 ^)/Fi K (m 2 ^) = (m 2 B - m 2 ^)/m 2 B for form factors. 
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Since asymptotically $^(77) =677(1 — 77), the logarithmic divergence of the 77 integral in Eq. 
( P-5| ) implies that the spectator interaction is dominated by soft gluon exchanges between 
the spectator quark and the strange or anti-strange quark of <fi. Hence, QCD factorization 
breaks down at twist-3 order. Note that the hard gluon exchange in the spectator diagrams 
is not as hard as in the vertex diagrams. Since the virtual gluon's momentum squared there 
is k 2 = (—ppB + VPk) 2 ~ —pfjm 2 B ~ Ph m b, where ph is the hadronic scale ~ 500 MeV, we 
will set a s ~ a s (y/phjnb) in the spectator diagrams. For the second term in (|3T5|), due to 
the end point divergence, the scale may correspond to a softer scale p s . However since a s p x 
is weakly scale-dependent, we can treat it at the ^Jp^m^ scale. The corresponding Wilson 
coefficients in the spectator diagrams are also evaluated at the ph scale. 

The infrared divergence is manifested in the integral Jq 1 df\lf\. However, it is known that 
the collinear expansion cannot be correct in the end point region owing to the transverse 
momentum (kr) of the quark which is averagely about 300 MeV, the order of the meson's 
size. Thus the lower limit of Jq 1 ^77/77 should be approximately proportional to 2{kj) /^ib] or 
equivalently, Jq dfj/fj can be approximately replaced by Jo dfj/(fj + (2k T )/m b ). A consistent 
treatment of kx in the calculation is still an issue since fcr itself is a higher twist effect in 
the QCD factorization approach. Thus we will treat the divergent integral as an unknown 
"model" parameter and write 

K./'? = l n (!+,„), (3.9) 

Jo 77 \Ph ) 

with ph being a complex number whose phase may be caused by soft rescattering [14|. 



We see that although the scattering kernel induced by the twist-3 LCDA of the kaon is 
formally power suppressed in the heavy quark limit, it is chirally enhanced by a factor of 
(2/i x /Aq CD ) ~ O(10), logarithmically enhanced by the infrared logarithms. 

Finally, we wish to remark that the leading-twist LCDAs of the <ft meson are given by 



(0(p,a)|^) 7mS (o)|o) = /^- F _p M y de e ^-$f(o, 



(<P(P,X)\s(x)a, u s(0)\0) = -ifl{ef^P v -el^P,) £ d£<* p -Qi{£), (3.10) 

where e* is the polarization vector of 0, £ is the light-cone momentum fraction of the strange 
quark in 0, and /J are vector and tensor decay constants, respectively, but the latter is 
scale dependent. Although $j| and <3>^ have the same asymptotic form, it is found that the 
transverse DA does not contribute to // and fu if light quarks are massless. The contribution 
of $^ to vertex corrections is suppressed by a factor of m^/rriB and hence can be neglected. 
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IV. ANNIHILATION AMPLITUDES 



As shown in the annihilation amplitude is formally power suppressed by order 
Nevertheless, it has been stressed in the PQCD approach that annihilation con- 
tributions in hadronic charmless B decays are not negligible ||16|| . There are four weak 
annihilation diagrams as depicted in Fig. 2. We first consider the annihilation amplitudes 
induced by (V — A)(V — A) operators. The first two diagrams, Figs. 2a and 2b, are fac- 
torizable diagrams and their contributions are of order m^/m 2 B and hence can be neglected. 
Indeed, the factorizable annihilation amplitude should vanish in — > limit owing to cur- 
rent conservation. It is easily seen that only o dd operators contribute to the nonfactorizable 
annihilation diagrams Figs. 2c and 2d. It turns out that the nonfactorizable annihilations 
are dominated by Fig. 2d owing to an endpoint contribution. Explicit calculations yield [see 
Eq. Q]: 



-2H< 



dpdidr] — — -™ + / dpdidr]— =- r - 



2H{Q(Y' - 1) 



V 



v(€v -€-v) 



where 



= as_ Cf_ 4tt 2 /b/k 



4tt jV fi iV, 



m 



le 'Pi 



(4.1) 



(4.2) 



and we have applied the approximation p « 1 and p = 1 — p ~ 0. In Eq. ( J4. 1|) the first term 
in brackets comes from Fig. 2d, while the second term from Fig. 2c. Since the soft phase of 
annihilation diagrams is not necessarily the same as that of the spectator diagram, we write 



Y' 



1 drj 
o fj 



In 



m B 



(4.3) 



where the phase is characterized by the complex parameter 

For (V — A)(V + A) operators O5 — Og, the twist-2 kaon DA makes no contribution. 
Therefore, we need to consider the twist-3 kaon LCDAs and $^ with the former being 
defined in the pseudoscalar matrix element ||13|| : 



(K-(P)\s(0)i l5 u(x)\0) 



fKfn 2 K 



m s + m u Jo 

The factorizable annihilation amplitude has the expression 



dfj^^ifj). 



(4.4) 



A f = AN C H 



2AN r H 



2 Vx 
m B 



Jo 




4^ 2 



Y '(y' 



1. 

2' 



(4.5) 
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(a) 



(b) 





FIG. 2. Annihilation diagrams for B — > decays. 

where we have applied the LCDAs: = 6£(1 - £), $^(77) = 1 and $f (77) = 67/(1 - 77). 

Likewise, the nonfactorizable annihilations induced by the penguin operators O5 and O7 read 

f^^x\ ) aVfV -\\ 1 [ l AC f 1 A^^(C\fcK(^\ 2-7] [ 



^l 2 



(5) 6r < r - + z w ft (fa if-,, }- 



where the dominated first term in brackets stems from Fig. 2d. Note that we have introduced 
a cutoff to regulate the infrared divergence occurred in the second term. 

Although the annihilation amplitudes Aj and A\* are formally of order (Aqcd/^) 2 , they 
receive two large enhancements: one from the chiral enhancement (2/i x /Aqcd) ~ 0(10) 
and the other from the logarithmic endpoint divergence of the infrared divergent integral 
Y' . Consequently, the annihilation effects can be sizable. Physically, this is because the 
penguin-induced annihilation contributions are not subject to helicity suppression. 

V. RESULTS AND DISCUSSIONS 



To proceed for numerical calculations, we employ the meson LCDAs as follows: 

/i 2 ) = 6r7(l - 77) (l + £ a^(/i 2 )C 2 3 / 2 (2r7 - 1)) , 



n=l 



$f(p) = A^p 2 (l-p) 2 exp 
d^(£) = $jf (£) = 6£(1 - £) 



1 / pm B \ 2 

2 V oj b J 



(5.1) 



where C^J 2 are Gegenbauer polynomials and the values of the Gegenbauer moments are 



available in [13], u B = 0.25 GeV, and N B is a normalization constant. We use the decay 
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constants fx = 0.16 GeV, /# = 0.19 GeV, = 0.237 GeV, and the running quark masses: 
m b (m b ) = 4.40 GeV, m s (m b ) = 90 MeV m d (m b ) = 4.6 MeV and m u (m b ) = 2.3 MeV. The 
next-to-leading-order Wilson coefficients in NDR and HV 75-schemes are taken from 
Table XXII of 0; they are evaluated at p = m b (m b ) = 4.40 GeV and A§L = 225 MeV. For 



MS 

form factors we use F^ K (mfy = 0.38 as a benchmarked value. Note that F 1 s - fs '(m|) = 0.407 
in the Bauer-Stech-Wirbel model 0, while it is 0.37 in a QCD sum rule calculation fl~8f . 

For the parameters pu in ( |3.9| ) and pa in ( |4.3| ), in principle they may be complex due to 
final-state soft rescattering. We find that the decay rate is much more sensitive to pa than 
to pu- Presumably, some information on the parameter p can be extracted from the study of 
B — > Kir modes. It has been shown recently in Jl4| that increasing the parameter \p&\ from 



1 to 2 would increase the corresponding error on the Kit branching ratios in which case it 
would require considerable fine-tuning of the strong interaction phase of Y' in annihilation 
diagrams to reproduce the experimental value of the branching ratio. Hence, it is reasonable 
to assume that the model parameters are in the range \p\ < 1. Writing p& = \pa\ exp(i8), 
the branching ratio of B — > <pK vs. the phase 5 is plotted in Fig. 3. We obtain 

B(B~ -> <j>K~) = (4.3±f;2) x 10" 6 , B(B° -> <f>K°) = (4.0±?;|) x 10" 6 , (5.2) 

where the central value corresponds to the default values Pa = Ph = and the errors 
come from the variation of \ph\ and \pa\ from to 1; that is, the theoretical uncertainties 
come from power corrections of twsit-3 spectator interactions and annihilation contributions. 
Therefore, the predicted branching ratio is consistent with CLEO and BaBar numbers, but 
smaller than the BELLE result [see ( |1.2| ) and ( |1.3| )1- The corresponding absolute ratio of 
the annihilation to penguin amplitudes depends on the annihilation phase and is at most of 
order 0.25. In the absence of annihilation effects, the branching ratios are given by 

B(B± - 4>K±) = ( ^gf^ ) ( 3 - 8 ± 0-6) x 10- 6 , 

/ pBK ( m 2 \ \ 

B(B° -> 0K°) = \ 00 * (3.6 ± 0.6) x 10" 6 , (5.3) 
V 0.38 J 

where the error arises from the variation of pn from to 1. 

Needless to say, the major theoretical uncertainty stems from the unknown model pa- 
rameters pu and pa- It should be stressed that the infrared divergence here is always of 
the logarithmic type and other possible linear divergence occurred in annihilation diagrams 
with twist-3 wave functions are explicitly canceled out. As stressed in passing, the infrared 
divergence stems from the misuse of the collinear expansion in the end point region where 
the effect of the quark's transverse momentum is important. Since kj< is a higher-twist effect 
in QCD factorization, at present we treat the infrared divergent integral in a model manner. 

Several remarks are in order, (i) The calculations are rather insensitive to the unitarity 
angle 7 as the CKM matrix element V ub V* s is considerably suppressed, (ii) The scattering 
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FIG. 3. Branching ratio of B~ — > <pK~ vs. the phase of the complex parameter pa [see Eq. 
where the dark (horizontal) and light bands correspond to \pa\ = ; 1> respectively, with 
the variation of \pn\ from to 1. 



kernel fjj is dominated by the twist-3 effect. However, since the Wilson coefficients Cn and 
C2«_i have opposite signs, it turns out that the magnitudes of 03-10 [see Eq. ( |3.2| )1 are slightly 
reduced by the twist-3 terms and therefore the branching ratio is suppressed by about 10% 
in the presence of twist-3 effects in spectator interactions, (iii) In the QCD factorization 
approach, the strong phase of the annihilation amplitude is of order a 2 s since it comes from 
the annihilation diagrams in which the gluon line is inset with an enclosed quark loop, 
resembling the vacuum-polarization bubble. Consequently, the phase of the annihilation 
contribution is likely dominated by the soft one induced by soft scattering as characterized 
by the parameter pa- This is in contrast to the PQCD approach where the annihilation 



contributions have large strong phases |16] 



VI. CONCLUSIONS 

We have analyzed the decay B — > <f)K within the framework of QCD-improved factoriza- 
tion and taken into account some of power-suppressed corrections. Our conclusions are: 

1. Although the twist-3 kaon distribution amplitude dominates the spectator interactions, 
it will suppress the decay rates of B — > (f)K slightly by about 10%. In the absence 
of annihilation contributions, the branching ratio is (3.8 ± 0.6) x 10 -6 for <pK~ and 
(3.6 ± 0.6) x 10~ 6 for (pK°. 
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2. The weak annihilation diagrams induced by (S — P)(S + P) penguin operators, which 
are formally power-suppressed by order (Aqcd / m b) 2 , are chirally and logarithmically 
enhanced. Therefore, these annihilation contributions are not subject to helicity sup- 
pression and in principle can be sizable. 

3. The branching ratio is predicted to be (4.3l?;2) x !0~ 6 for B~ -> (pK~ and (4.0lf;|) x 
1CT 6 for B° — > (j)K°, where theoretical uncertainties come from power corrections 
of twist-3 spectator interactions and annihilation contributions. The corresponding 
absolute ratio of annihilation to penguin amplitudes depends on the annihilation phase 
and is at most 25%. 
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